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In the last few years the Fermi-LAT instrument has detected GeV gamma-ray emission from a few 
novae. Such GeV emission can be interpreted in terms of an inverse Compton process of electrons 
accelerated in a shock. It is expected that hadrons can be accelerated in the same conditions, but 
reaching much higher energies. They can produce a second component in the gamma-ray spectrum at 
TeV energies. We performed follow-up observations of selected novae and dwarf novae in a search of 
the second component in the gamma-ray spectrum. This can shed light on the acceleration process 
of leptons and hadrons in nova explosions. We have performed observations with the MAGIG 
telescopes of 3 sources, a symbiotic nova YY Her, a dwarf nova ASASSN-13ax and a classical nova 
V339 Del shortly after their outbursts. 


I. INTRODUCTION 

A classical nova is a thermonuclear runaway leading 
to the explosive ejection of the envelope accreted onto 
a white dwarf (WD) in a binary system in which the 
companion is either filling or nearly filling its Roche 
surface [H-Q . They are a type of cataclysmic variables, 
i.e. optically variable binary systems with a mass 
transfer from a companion star to WD. Novae are 
typically detected first in optical observations when 
the brightness of the object increases by 7-16 magni¬ 
tudes. The energy spectra of novae often contain a 
thermal X-ray continuum. The symbiotic novae, like 
the classical novae, are also initiated by a thermonu¬ 
clear explosion on the surface of the WD. However 
in the case of symbiotic novae, the WD is deep im¬ 
mersed in the wind of a late-type companion star (see 
e.g. [1,[1). 

The diffusive shock acceleration at the blast wave 
of symbiotic novae was expected to accelerate parti¬ 
cles up to energies of a few TeVs @. In 2010 the 
first GeV gamma-ray emission was detected by the 
Fermi-LAT from the symbiotic nova V407 Cyg Q- 


Such gamma-ray emission can be explained in terms 
of either leptonic or hadronic models 0 , [ 1 ]. In the 
former case, the local radiation fields create a target 
for the inverse Compton (IC) scattering of the elec¬ 
trons. On the other hand, protons accelerated in the 
same conditions can interact with the matter produc¬ 
ing gamma-rays via proton-proton interactions. Sev¬ 
eral models have been put forward to explain the GeV 
radiation. For instance, the GeV gamma-ray emission 
can be attributed to the IC process on the strong ra¬ 
diation field of the red giant 0. The same model 
predicts a second component in the TeV range due to 
proton-proton interactions with the wind of the red 
giant. Also [l^ consider acceleration of leptons and 
hadrons in the nova shock. In that model the magnetic 
field, which determines the acceleration efficiency, is 
obtained assuming an equipartition with the thermal 
energy density upstream of the shock. The GeV y-ray 
emission is then a product of IC scattering of the nova 
light by the electrons. 

In the last few years the Fermi-LAT has discov¬ 
ered GeV gamma-ray emission from a few more no¬ 
vae: V1324 Sco, V959 Mon, V339 Del, and V1369 
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Cen dini. Most of these sources are classical novae. 
Contrary to the symbiotic ones, they do not exhibit a 
strong wind of the companion star. Interestingly, sym¬ 
biotic and classical novae all exhibit similar spectral 
properties of the GeV emission. In classical novae the 
particles acceleration can occur e.g. on a bow shock 
between the nova ejecta and the interstellar medium 
or in weaker internal shocks due to inhomogenuity of 
the nova ejecta Q. In particular orbital motion of the 
system can lead to shaping the nova ejecta into a com¬ 
bination of a faster polar wind of the WD ejecta, and 
a denser material drifted along the equatorial plane 


So far no very-high-energy (VHE; E>100 GeV) 
gamma-ray emission has been detected from any nova 
event. VERITAS observations of V407 Gyg which 
started 10 days after the nova explosion yielded a 
differential upper limit on the flux at 1.6 TeV of 
2.3 X 10~^^ ergcm~^ s~^[T^ 

Beginning in Eall 2012 the MAGIG telescopes con¬ 
ducted a nova follow-up program in order to detect a 
possible VHE gamma-ray component. The program 
was first aimed on symbiotic novae. After the reports 
of detection of GeV emission from classical novae by 
the Eermi-LAT, the program was extended also to 
bright classical novae and opened up to additional out¬ 
bursts from cataclysmic variables. 

In here we report on MAGIC and Fermi-LAT 
(see Section |TT] for description of the used instru¬ 
ments and analysis methods) observations of the 3 
sources observed within this program: a symbiotic 
nova YY Her f Section IHI Al) . a dwarf nova ASASSN- 
13ax (Section IIII HI) and a classical nova V339 Del 
ISection IHI Cl) . 


II. INSTRUMENTS 


data were analyzed using the standard analysis chain 
[isl . . The significance of a gamma-ray excess was 
computed according to Eq. 17 of El. The upper lim¬ 
its on the flux were calculated following the approach 
of using 95% C.L. and accounting for a possible 
30% systematic uncertainty on the effective area of 
the instrument. 


B. Eermi-LAT 

The Fermi-LAT, launched in 2008, is a pair- 
conversion telescope that detects photons with en¬ 
ergies from 20MeV to > 300 GeV E^- Thanks 
to a large field of view (^ 2.4 sr), the Fermi- 
LAT observatory, operated in scanning mode, pro¬ 
vides coverage of the full sky every three hours en¬ 
abling searches for transient sources and overlap with 
ground-based observatories. We analyzed the LAT 
data in the energy range 100 MeV — 300 GeV us¬ 
ing an unbinned maximum likelihood method 
as implemented in the Fermi Science Tools v9r32p5, 
the P7REP_SOURCE_V15 LAT Instrument Response 
Functions (IRFs), and associated standard Galactic 
and isotropic diffuse emission models [l^. We selected 
events within a region of interest (ROI) of 15° cen¬ 
tered on the LAT best position (following Q) for V339 
Del and required a maximum zenith angle of 100° in 
order to avoid contamination from Earth limb pho¬ 
tons. Additionally, we applied a gtmktime filter (no.3) 
recommended for combined survey and pointed mode 
observations [i^, selecting good quality data at times 
when either the rocking angle was less than 52° or the 
edge of the analysis region did not exceed the max¬ 
imum zenith angle at 100°. Sources from the 2FGL 
catalogue [2l| located within the ROI were included 
in the model used to perform the fitting procedure. 


The three sources were first detected and observed 
by optical instruments. The results of the MAGIC HI. SOURCES OBSERVED 

observations were supported by the analysis of quasi- 

simultaneous Fermi-LAT observations. We report here results of the MAGIC and Fermi- 

LAT observations of YY Her, ASASSN-13ax and V339 
Del. 

A. MAGIC telescopes 


The VHE gamma-ray observations were performed 
with the MAGIC telescopes. MAGIC is a system of 
two 17 m Cherenkov telescopes located on the Canary 
Island of La Palma at a height of 2200 m a.s.l. The 
telescopes can perform observations of gamma rays 
with energies as low as ~50GeV and up to tens of 
TeV. During Summer 2011 and 2012 MAGIC under¬ 
went a major upgrade E3|- After the upgrade the 
sensitivity of the MAGIC telescopes in the best en¬ 
ergy range (> 300 GeV) is 0.6% of Crab Nebula 
flux in 50 h of observations E3- All the data used 
for this paper were taken after the upgrade. The 


A. YY Her 

YY Her is a symbiotic nova system that undergoes 
a recurrent pattern of outbursts. MAGIC observa¬ 
tions of YY Her occurred on the night of 2013 Apr 
22''‘^/23’'‘^, 7 days after the optical maximum. No sig¬ 
nificant VHE gamma-ray emission was detected. We 
computed flux upper limits at 95% confidence level ob¬ 
taining < 5.0x 10“^^phcm“^ s“^ above 300 GeV. Also 
in Fermz-LAT no emission was detected over a longer 
interval 2013 Apr lO*^' to Apr 30*** (MJD 56392.5 to 
56412.5). Upper limits at 95% confidence level were 
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FIG. 1: Differential upper limits on the flux from YY 
Her as measured by the Fermi-LAT (empty squares) and 
MAGIC (full circles). See text for details of the time 
ranges covered by the points. For comparison a spectrum 
of Crab Nebula is shown with a gray curve. 


FIG. 2: Differential upper limits on the flux from 
ASASSN-13ax as measured by the Fermi-hNT (empty 
squares) and MAGIC (full circles). See text for details 
of the time ranges covered by the points. For comparison 
a spectrum of Crab Nebula is shown with a gray curve. 


set as 2.8 x 10“®phcm“^ s“^ above 100 MeV. Differ¬ 
ential upper limits obtained from the Fermi-LAT and 
MAGIC observations of YY Her are shown in Fig. [T] 


B. ASASSN-13ax 

ASASSN-I3ax is a member of a different class of 
cataclysmic variables, the dwarf novae, which are 
known for significantly weaker optical outbursts (2-6 
magnitudes) than classical novae. Instead of under¬ 
going a thermonuclear explosion on the surface of the 
WD, these outbursts are caused by the gravitational 
energy release from a partial collapse of the accre¬ 
tion disk surrounding the WD. The MAGIC observa¬ 
tions were performed on two consecutive nights start¬ 
ing on 2013 Jul 4**^, soon after the optical outburst 
seen on 2013 Jul 1®*. In the absence of detectable 
VHE emission, upper limits at 95% confidence level 
were set as < 1.5 x 10“^^phcm“^ s“^ above 300 GeV. 
Emission was not detected in the LAT over the in¬ 
terval 2013 Jun 25“^ to Jul 15*^ (MJD 56468.5 to 
56488.5). Upper limits at 95% confidence level were 
set as 1.6 x 10“®phcm“^ s“^ above 100 MeV. Differ¬ 
ential upper limits obtained from the Fermi-LAT and 
MAGIC observations of ASASSN-13ax are shown in 
Fig.H 


C. V339 Del 

V339 Del was a fast, classical nova detected by op¬ 
tical observations on 2013 Aug 16*^ (CBET #3628). 
The nova was exceptionally bright reaching a magni¬ 
tude of 5 mag (see top panel of Fig. [3]), and it 


triggered follow-up observations at frequencies rang¬ 
ing from radio to VHE gamma-rays. Photometric 
measurements suggest a distance for V339 Del of 
4.5±0.6kpc [l^l. The spectroscopic observations per¬ 
formed on MJD 56522.1 revealed emission wings ex¬ 
tending to about ±2000 km/s and a Balmer absorp¬ 
tion component at a velocity of 600 ± 50km/s |^ . 
The pre-outburst optical images revealed the progen¬ 
itor of nova V339 Del to be a blue star [^ . 

Originally MAGIC observations of V339 Del were 
motivated by its extreme optical outburst. Soon af¬ 
ter MAGIC started observations they were addition¬ 
ally supported by the detection of GeV emission by 
the Fermi-LAT from the direction of V339 Del. The 
MAGIC observations started already on the night of 
2013 Aug 16/17*^*', however they were marred by bad 
weather conditions. The good quality data used for 
most of the analysis spanned 8 nights between 2013 
Aug 25**^ and Sep 4*^*^. The total effective time was 
11.6 h. In addition to the nightly upper limits we per¬ 
formed a dedicated analysis of the poor quality (af¬ 
fected by calima, a dust layer originating from Sa¬ 
hara) night of 2013 Aug 16/17*^. We applied an es¬ 
timated energy and collection area corrections based 
on LIDAR measurements [1^. No VHE gamma-ray 
signal was found from the direction of V339 Del. We 
computed a night by night integral upper limit above 
300 GeV (see bottom panel of Fig. The differen¬ 
tial upper limits for the whole good quality data set 
computed in bins of energy are shown in Fig. 01 

Nova V339 Del was the subject of a Fermi Tar¬ 
get of Opportunity (ToO) observation [2^ triggered 
by the optical discovery (CBET #3628); the ToO 
started on 2013 Aug 16*^ and lasted for six days. The 
gamma-ray emission from V339 Del was first detected 
by Fermi-hAT in 1-day bins on Aug 18*^ Q. The 
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FIG. 3: Multi-wavelength light curve of V339 Del during 
the outburst in August 2013. Top panel: optical observa¬ 
tions in the V band obtained from AAVSO-LCGQ service. 
Middle panel: the Fermi-hAT flux (filled symbols) and up¬ 
per limits (empty symbols) above 100 MeV in 1-day (cir¬ 
cles, thin lines) or 3-day (squares, thick lines). A 95% C.L. 
flux upper limit is shown for time bins with TS<4. Bottom 
panel: Upper limit on the flux above 300 GeV observed 
with MAGIG telescopes. The gray band shows the obser¬ 
vation nights with MAGIG. The dashed gray line shows a 
MAGIC observation night affected by bad weather. 

“http://w ww.aavso.org/leg 



FIG. 4: Differential upper limits on the flux from V339 
Del as measured by MAGIC (filled circles) and the flux 
measured by Fermi-hAT (empty crosses) in the same time 
period, 25*** of August and 4*^ of September. For compari¬ 
son a spectrum of Crab Nebula is shown with a gray curve. 


bins with TS<9. In the period coincident with the 
MAGIC observations (2013 Aug 25*^ to Sep 4*^) the 
Fermi-hAT spectrum can be described by an effec¬ 
tive power law with an index of 2.37 ± 0.17 and flux 
above 100 MeV of (0.15 ± 0.04) x 10“®phcm“^s“^. 
The rather low statistical significance (TS=49) does 
not constrain the value of an exponential cut-off of the 
emission in this period. Note, however, that the most 
energetic photon, with E = 5.9 GeV was recorded 
on Aug 30**^, i.e. within the time period covered by 
the MAGIC observations. The Fermz-LAT analysis 
for a broader time range, 2013 Aug 22"'^ to Sep 12**' 
(MJD 56526-56547), covering the whole decay phase 
of the Fermi-LAT light curve allowed us to obtain a 
more significant signal with a TS of 121. Neverthe¬ 
less we obtain a similar value of flux above 100 MeV, 
(0.13 ± 0.03) X 10“®phcm“^ s“*, for this broader pe¬ 
riod. The spectrum in this case can be described, with 
improved significance of 3.3 (t with respect to the sim¬ 
ple power law, by an exponentially cut-off power law 
with an index of 1.44 ± 0.29 and a cut-off energy of 
1.6 ±0.8 GeV. 


emission peaked on Aug 22'“* and entered a slow de¬ 
cay phase afterwards (see middle panel in Fig. [3]). For 
the light-curves, the data were fit using a power law 
spectral model initially leaving the photon index and 
the normalization free to vary. We then fixed the pho¬ 
ton index at the average value of 2.3 calculated over 
the most significant detections (Test Statistic values 
TS>9)[2^. The LAT Spectral Energy Distribution 
(SED) of V339 Del shown in Fig. 0] was extracted 
in five logarithmically spaced energy bins from 100 
MeV to 100 GeV. Similarly to the light-curves, energy 
binned data shown in Fig. 0] were fit using a simple 
power law and showing a 95% G.L. upper limit for 


IV. CONCLUSIONS 

The MAGIC telescopes performed observations of 3 
objects: the symbiotic nova YY Her, the dwarf nova 
ASASSN-13ax and the classical nova V339 Del. No 
significant VHE gamma-ray emission was found from 
the direction of any of them. Out of these three ob¬ 
jects, V339 Del is the only one detected at GeV ener¬ 
gies. It has also extensive optical observations which 
shed some light on both the companion star and the 
photosphere of the nova. Therefore it has the highest 
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potential for constraining the leptonic and hadronic 
processes in novae. MAGIC will continue follow-up 
observations of the promising novae candidates in the 
following years. 
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